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ABSTRACT 

An exact, axially symmetric solution to the Einstein-Klcin-Gordon field equations is 
employed to model the dark matter in spiral galaxies. The extended rotation curves from a 
previous analysis are used to fit the model and a very good agreement is found. It is argued that, 
although our model possesses three parameters to be fitted, it is better than the non-relativistic 
alternatives in the sense that it is not of a phcnomcnological nature, since the dark matter 
would consist entirely of a scalar field. 
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Introduction. 



Since the pioneering works by Oort and Zwicky, back in the 1930's ( Port 193^ ; Zwicky 1933), the 



existence of dark matter in the Universe has been firmly established by astronomical observations at very 
different length-scales, ranging from the neighbourhood of the Solar System to the clusters of galaxies. 
What it means is that a large fraction of the mass needed to produce, within the framework of Newtonian 
mechanics, the observed dynamical effects in all these very different systems, is not seen. This puzzle has 
stimulated the exploration of lots of proposals, and very imaginative explanations have been put forward, 
from exotic matter to nonrelativistic modifications of Newtonian dynamics. 



In particular, the measurement of rotation curves (RC) in galaxies shows that the coplanar orbital 
motion of gas in the outer parts of these galaxies keeps a more or less constant velocity up to several 
luminous radii. The discrepancy arises when one applies the usual Newtonian dynamics to the observed 
luminous matter and gas, since then the circular velocity should decrease as we move outwards. The most 
widely accepted explanation is that of a spherical halo of dark matter, its nature being unknown, which 
would surround the galaxy and account for the missing mass needed to produce the flat RC. Another 
possibility, considered much less often, is the so called Modified Newtonian Dynamics (MOND), which 
was put forward by Milgrom (1983); in this model, the usual second Newton law would broke at "small" 
accelerations, as compared to some (in principle universal) acceleration parameter, a . Although it seems to 
provide for a very good phenomenological description of the RC, it lacks, at least until now, a more sound 
theoretical basis. 



Our aim here is to give another explanation to the dark matter problem in spiral galaxies, this time 
using a fully relativistic approach, and making use of the well known scalar fields. Combining an exact 
solution of the Einstein-Klcin-Gordon field equations with the observations from luminous matter and gas, 
we are able to reproduce the flat extended RC of spiral galaxies. Scalar fields are the simplest generalization 
of General Relativity and they can be introduced on very fundamental grounds, as in the case of the 
Brans-Dicke, the Kaluza-Klein and the Super-Strings theories. They appear also in cosmological models, 



like inflation, and in general in all modern unifying theories. Recently, it has been suggested (Cho & Keum 
1998| ; [Peebles 1999|) that a massive scalar field could account for the dark matter at cosmological scales, and 



a previous work ( |Guzman fc Matos 200C ) has shown a preliminary analysis in the context of spiral galaxies. 
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This is encouraging since it makes models like the one proposed here to seem more plausible. 

The paper is organized as follows: in we introduce the field equations and the explicit solution for 
an axially symmetric configuration; in the geodesic equations are written and the model for the dark 
matter in a spiral galaxy is introduced; section §^ gives the main results concerning the fitting of the model 
to the observations, and in §|| some concluding remarks are done. Finally, two appendices describe some 
geometrical properties of our metric and make a brief conceptual comparison with the dark halo and the 
MOND hypotheses. 



2. The field equations and their solution. 

As mentioned in the Introduction, scalar fields appear in a natural way within the framework of 
unifying theories. As example we mention Kaluza-Klein (KK) and Super-Strings (SS), where the scalar field 
appears in the effective action after dimensional reduction. Let us begin with the most general scalar-tensor 
theory of gravity, as given by the action: 



S = J d 4 aV-3 (^-■^ J R+ (V$) 2 G + f($)J . (1) 
where F, V and G are functions of only, g is the determinant of the metric g^ and k q = !6irG/c . We 



can perform a conformal transformation to some other frame by means of the redefinitions (Damour & 



Esposito-Farese 1992| ; [Frolov 1999|) 



(2) 
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(3) 



where the prime denotes d/d&. The action then takes the form: 



S = [ d*xV=d[-- + 2(V0) 2 + vm, 
J K 



(4) 
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where R is the four dimensional scalar curvature and g the determinant of the metric g^ v . <f> and V((j>) 
repectively are the scalar field and the scalar potential in the new frame, and all the other quantities are 
also calculated using the new metric. The choice of the Einstein frame, implicit in this form of the action, 
is made because the field equations are more easily solved in this frame. As is well known, since the 
coupling of a scalar field with gravity is defined up to a conformal transformation, there is some ambiguity 



in regarding a specific frame as the "physical" frame (see, e.g., Damour & Esposito-Farese 1993). We shall 
consider in detail what happens both to the action above and to the equations of motion in section 0. 



The next step is to decide what kind of potential V(<f>) is the most convenient for modeling a galaxy. 
Let us reason as follows. It is known that the energy density of the dark matter in the halo of the galaxies 
goes like 1 /r 2 . The energy momentum tensor of the scalar field is basically the sum of quadratic terms of 
the scalar field derivatives plus the scalar potential, i.e., ~ 4>,n4>.u + V(<f>) ~ 1/r 2 . If we assume that the 
term (f> a<f> v ~ 1/r 2 as well as the term V(<f>) ~ 1/r 2 , from the first assumption we infer that <f> ~ l n (r), and 
from the second assumption we arrive at V(<f>) ~ exp(— 2<j>). So, in what follows we will take the potential 
V{(j>) = Aexp(— 2a(f>) and we shall consider the four dimensional action: 



S= [d 4 x^g~[- — + 2(V0) 2 +e- 2a *A]. 
J ^0 



For the time being, from this very general setting we shall look for an exact solution to the field 
equations that will serve us as a model for a spiral galaxy. Now, since the velocity of the gas and the 
red shift measurements in a galaxy are made over the equatorial plane, it is reasonable to impose axial 
symmetry on the solution we are looking for, in contrast with the usual spherical dark halo profile used 



in most studies ( Begeman et al. 1991 ). Moreover, the fact that a substantial amount of the total mass in 
these galaxies is in the form of dark matter suggests that, in a first approximation, the observed baryonic 
mass (both stars and gas) will not contribute significantly to the total energy density of the system, at least 
in the region outside the luminous disk; instead, the scalar matter will determine the space-time curvature, 
and the material particles will move on geodesies determined (almost) by the energy density of the scalar 
field. Finally, the RC exhibits a constant velocity of the order of 100-350 km/s, which, compared to the 
velocity of light, clearly allows one to consider the galaxy as a static system. 



With the above simplifying assumptions, the most general axially symmetric, static metric can be 



written as: 



ds 2 = j[e 2k (dzdz) + W 2 d<p 2 } - fc 2 dt 2 , (5) 

where z := p + i£, the bar means complex conjugation, and the real valued functions /, W and k depend 
only on z and z (or equivalently on p and C) . 

After varying the action in equation ([!]), one obtains the following field equations: 

i?^ = Ko [2<f>,ii<f>,v + ^g^V {(/))}, (6) 

which are the Klein-Gordon and Einstein field equations, respectively, and (i, v = 0, 1, 2, 3. A very powerful 
technique, known as the harmonic maps ansatz, can be employed to find families of solutions to the 
equations @, starting with the metric (||). The details can be found in Matos (1989; 1994; 1995) and 
Guzman & Matos (1999), so we shall only describe it very briefly here. 



2.1. The harmonic maps ansatz. 

In a few words, the main idea behind the method is to re-parameterize the functions in metric (^|) with 
convenient auxiliary functions which will obey a generalization of the Laplace equation, along with some 
consistency relationships; the latter are usually quite difficult to fulfill, and great care and intuition must 
be taken in order to get a system of equations both workable with and interesting enough. In this case we 
shall take / = e A , and assume that / and <j> are functions of W, which in turn is a function of z and ~z alone. 
After lengthy but straightforward calculations one is left with the system: 
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W W n 

k z = ^- z +—\ 2 z W^ + n a W4, 2 z W, Zl (7) 

and a similar equation for k >z , replacing z by z. The symbol A stands for a generalized Laplace operator, 
such that for every function h(z,~z), Ah :— (Wh tZ ),z + {Wh t g) >z . 



2.2. The model for the galaxy. 

Regarding the set of equations (^|), it can be noted that the last equation and its complex conjugate 
are integrable once the functions A, <fr and W have been integrated. The first three equations, however, are 
highly coupled since \J—g — W^exp (2k — A)/2; moreover, the operator A itself contains W . In spite of this, 
we have been able to obtain a not too restrictive solution, which can be written as: 



A = InM + ln/o, 
1 



V 



4/o 



TnM, 



K M' 

e 2fc = M iZj M, (8) 
where /o and 4>o are integration constants with 

e -2^o = 4/ A/ Ko (9) 
and M = W is, as stated before, a function of z and z, restricted only by the condition 

MM,- Z = M tZ Mj, (10) 

but otherwise arbitrary. Observe that A and kq are fundamental constants of the theory but fo and 4>o are 
integration constants, i.e., they are different for each space-time (each galaxy), fulfilling the relation (|^). 
The values of these constants will determine the characteristics of a particular galaxy. 
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We shall take the solution (||) as the general relativistic description of the galaxy, considering a 
particular choice of M, namely, M = zz/rQ, where ro is a constant with dimensions of length. 



Since the particles compossing the gas from which observations arise are small compared to the whole 
galaxy, they can be considered as test particles moving on the background metric (Q). Therefore, the next 
step is to study the geodesic motion of test particles along the equatorial plane. From metric (^|) we can 
write, for material particles: 




As the solution is axially symmetric and static, there will be two constants of motion, namely, the angular 
momentum per unit mass, 



where r is the proper time of the test particle. In order to obtain useful information from these constants, 
it is convenient to write the line element as: 



3. Geodesic motion along the equatorial plane. 



_ W 2 dip 



and the total energy of the test particle, 





(12) 



(13) 



since the squared three- velocity, v 2 



is given by: 



■2 
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where a,b = 1,2,3. On the other hand, for a material freely falling observer (i.e., an observer in geodesic 
motion) we must have ds 2 = —c 2 dr 2 , and equating this expression with equation (|l2|), we can arrive to: 



c 4 f 2 

(14) 



/ — v 2 /c 2 

If we now identify the equatorial plane of the galaxy with the plane C = 0, the geodesic equation ( [LI 
reduces to the following, after inserting the constants of motion A and B: 



lo t (dp\ J , B 2 f A' _, 



r"\t) + wi-sf=^ <15 » 

This equation describes the motion of test particles in the equatorial plane of the galaxy, and in every 
particular trajectory the constants A and B remain so, i.e., constant. However, a key point here is that, 
shall we change of trajectory, the values for A and B will change accordingly for the new trajectory. Since 
the RC give an average of the circular velocity of particles in the galaxy, we shall consider circular orbits 
only, for which p = 0, so that v in equation ( |l2] ) can be identified with v c = v C i rcu i ar . From this, an 
expression for B in terms of v 2 can be written down: 



B 



and since, as stated before, v <§; c, this gives: 



W 



f-v 2 /c 2 f 



(16) 

Using now the form of / given by equations (^), one arrives at the following remarkably simple relation 
between v and B: 



VDM = f B, (17) 

where we have written vdm instead of v to stress the fact that this velocity for test particles is due to 
the scalar dark matter. Formula ( |l7j ) is the main result of this model, it states the way how the circular 
velocity due to the dark matter is determined by the angular momentum from each orbit. It is remarkable 
that formula ( [1^ ) is invariant under conformal transformation = Fg^ LV of the metric, i.e., this formula 
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is valid for the metrics g^, and g^ 



In order to gain some insight into the physical meaning of the solution (Q), we write it in Boyer- 
Lindquist-like coordinates (Schwarzschild coordinates) (r,6,ip), related to (p,(,tp) by p = \/r 2 + b 2 sin (9, 
C = rcosO; metric (||) then reads: 



., 1 + b 2 cos 2 6/r 2 ( dr 2 2, a 2\ r 2 + b 2 sin 2 8 2 2 r 2 + b 2 sin 2 9 , 



ds 2 = — f '— r—nn-2 + r 2 d9 2 + dip* - f c 2 : dt 2 . (18) 

foro \l + b 2 /r 2 J f r r 

On the other hand, the effective energy density Hum is given by: 



-V((j>) = — — — 

^ Ko(r 2 + b 2 sin 



Mbm = ^V{4>) = ., ,., . , ir (19) 



The fact that this energy density is negative does not constitute a serious drawback since, as mentioned 
before, we most perform a conformal transformation of the metric (|^) in order to obtain the action 
corresponding to a theory with a more physical interpretation. 



To be able to obtain more quantitative results and to compare this model with the most usual 
approaches, one further assumption must be made, regarding the constant of motion B. The observed 
luminous matter in a galaxy behaves in accord to Newtonian dynamics to a good approximation, so that 
its angular momentum per unit mass will be B = Vl X D, where vl is the contribution of the luminous 
matter, and D is the interval from the metric as written in equation (|T^); as we are in the equatorial 
plane, on circular orbits and at one particular instant, we have dp = d£ = dt = 0, and it is easy to check 
that D — y^(r 2 + b 2 )/ fo. It is now reasonable to substitute this value for B in equation (0), since the 
expression for vdm represents the velocity of test particles due to the presence of the scalar field; we get: 



vdm = hv L VW+Wh, (20) 

Noting that the total kinetic energy will be well approximated by the sum of the individual contributions, 
i.e., \J1mVQ ~ \j2mv 2 L + l/2mv 2 as + mv 2 DM , we arrive at the final form of the velocity along circular 
trajectories in the equatorial plane of the galaxy: 
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v 2 c=v 2 L (l + f (r 2 + b 2 )) + v 2 gas , (21) 

where the constants /o and 6 will be parameters to adjust to the observed RC. To this end we shall proceed 
as follows: we take the photometric and RC data for 6 spiral galaxies from Begeman et al. (1991) and 
Kent (1987), as listed in Table |. This data is fitted using a non-linear least squares routine adding a third 
parameter, namely, the usual mass-luminosity ratio M/L, which is taken to be constant in each particular 
galaxy; when there exist disk and bulge observations, two M/L ratios are assumed. The total luminous 
mass at a distance r from the center of the galaxy will be Ml(t-) = (M/L) x L(r), i.e.: 

vlir) = 2MiM (22) 



Combining equations (|2l| ) and ( |22| ), and including the 21 cm data from gas, we are led to: 



vUr) = GM^ {1 + Mr2 + b2)) + v2 ^ (23) 

In the next section |] we shall show how this compares with the actual observational data. It can also 
be useful to compare this derivation with more common explanations for the RC, two of which are briefly 
described in appendix H. 



4. Results. 

The main results are shown in Fig. ^ and in Table |[ Figure [| shows the observational RC (for 
simplicity we have omited the error bars) as well as the fitted curves using equation (|23|). Shown are also the 
individual contributions from luminous matter, gas and the scalar field. It can be noted that the agreement 
is quite good (within 5% in all cases), which could have been expected since there are three parameters 
to be adjusted. However, it should be noted that this approach is made on a very solid theoretical basis, 
because we have begun with a relativistic description of the galaxy. Moreover, the dark matter in this 
model would be entirely constituted by the scalar field. 



For the parameter M/L our results are in very good agreement with previous analyses which employ 
the dark-halo and the MOND approaches (e.g. Begeman et al. 1991 ). The remaining two parameters, /o 
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Fig. 1. — Curves of vc (km/s) vs r (kpc) for the observational data (dots) and the fits from equation (|23|) 
(solid). Also shown are the individual contributions from luminous matter (short dashed), the dilaton (long 
dashed) and gas (dotted). In the case of NGC 2841 and NGC 7331, the lower short dashed curve represents 
the contribution from the bulge of these galaxies. 
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and b, serve only to determine completely the metric, and they do not have a direct physical interpretation 
other than as part of the scalar field energy density. In Table || the best fit parameters are listed, along 
with the formal la error in the fitted parameter. 



In this work we have obtained an exact, axially symmetric and static solution to the field equations of 
gravity coupled with a scalar field. This solution has been successfully employed as a model for a spiral 
galaxy, and in particular, we have been able to reproduce the RC of matter in these galaxies with an 
excellent agreement, both with the observations themselves and with previous analyses of this kind of data. 

It should be stressed that, although our model has three parameters to be fitted, which in general 
allows for a great flexibility, they are obtained from a theory that is of a fundamental nature, namely, the 
low energy limit of a family of unification theories. This makes the calculations herein shown to be natural 
since no ad hoc hypotheses are needed, in contrast to the dark halo or the MOND models. We conclude 
that this work enables us to state that scalar fields are strong candidates to constitute the dark matter, not 
only at a cosmological scale, but also within spiral galaxies. 

We want to express our acknowledgment to the relativity group in Jena for its kind hospitality and 
partial support. This work was also partly supported by CONACyT, Mexico. 



In order to gain some insight into the solution (|18|), we shall consider a couple of issues related to the 
geometrical and topological aspects of the metric. By defining: 



5. 



Conclusions. 



A. Some geometrical aspects of the metric. 



Xi = 



In \/ r 2 + b 2 sin 2 9, 




x 3 = <p, 



x 4 = f t 



we have: 
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Table 1: Sample of galaxies 



Galaxy 


Type 


Distance 


Luminosity 


Rhi 






(Mpc) 


(1O 9 M ) 


(kpc) 


NGC 2403 


Sc(s)III 


3.25 


7.90 


19.49 


NGC 2903 


Sc(s)I-II 


6.40 


15.30 


24.18 


NGC 6503 


Sc(s)IL8 


5.94 


4.80 


22.22 


NGC 3198 


Sc(rs)I-I 


9.36 


9.00 


29.92 


NGC 2841 


Sb 


9.46 


20.50 


42.63 


NGC 7331 


Sb(rs)I-I 


14.90 


54.00 


36.72 



Table 2: Best-fit parameters 



Galaxy 


{M/L) dlsk 


{M/L) bulge 


b 

(kpc) 


fo 

(kpc- 1 ) 


NGC 2403 


1.75 




1.63 


0.0116 




0.04 




0.003 


5.7 x 10~ 4 


NGC 2903 


2.98 




8.33 


0.0043 




0.12 




0.03 


4.0 x 10~ 4 


NGC 6503 


2.12 




1.79 


0.013 




0.09 




0.01 


8.7 x 10~ 4 


NGC 3198 


2.69 




7.83 


0.0054 




0.08 




0.02 


3.0 x 10~ 4 


NGC 2841 


5.39 


3.25 


13.85 


0.0039 




0.34 


0.36 


0.16 


2.5 x 10~ 4 


NGC 7331 


5.06 


1.11 


0.845 


0.0013 




0.23 


0.06 


0.002 


8.9 x 10~ 5 
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ds 2 = j^re 2xi (dxj + dx\ + dxj - c 2 dx\). (Al) 

From this, it can be argued that our coordinates (9, ip) are not really angles, but in fact just cartesian 
coordinates. This, however, can only be stated if we know the global topology of our space-time, which 
we do not. As a counterexample, we can consider the 2-torus, T 2 ; in this case, the metric is not only 
conformally flat but in fact flat altogether, i.e., ds 2 — dx 2 + dy 2 ; with a proper rescaling, this can be written 
ds 2 = dd\ + d9\, where now 9\ and 02 vary from to 2ir and can be thought of effectively as angles. A 
similar example is the 2-sphere, S 2 . Therefore, our point of view is that, since interpreting (r 7 6,(p) as 
spherical-like coordinates allows us to reproduce the rotation curves for galaxies, we can consider them in 
such way, and in particular, ip does represent angles about the axial direction. 



B. Dark-halo profiles and the MOND hypothesis. 

The most commonly accepted approach to explain the RC is to assume that there is some kind of 
unseen matter around the visible part of the galaxy, forming what is usually called a 'halo'; then, a mass 
density profile for this dark matter is formulated and combined with data from the visible and 21 cm 
observations, and the model is fitted to the RC as obtained from red-shifts in the galaxy. A quite broad 
family of density profiles is given by ( |Zhao 1996| ): 



Phalo(r^ - 



(r/r e )i[l + (r/r c ) a ]W-l)/ a ' 

where po is the central density and r c is the 'core' radius, both of the halo. In particular, the simplest and 
most widely used profile is the so called modified isothermal sphere (MIS), for which (a,/3, 7) = (2,2,0). 
This model is attractive because there are cosmological arguments which seem to suggest that, under the 
suitable conditions, astronomical objects of this kind might actually evolve, the dark matter being cold 



or hot depending on the evolutionary arguments (see, e.g., Kravtsov et al. 1998 and references therein); 
however, a completely satisfactory evolution scenario remains to be derived. In the actual fitting procedure, 
it is usually more convenient to work with the asymptotic circular velocity Vh obtained from the isothermal 
sphere halo, by applying the virial theorem and Newton's law: 



v h = yJknGp r 2 c . 
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Beginning with the luminosity observations, L(r), there are three parameters to be adjusted (four in the 
case of galaxies with separate observations from the disk and the bulge): the ratio M/L, usually assumed 
to be constant over the whole optical disk, the core radius r c and the asymptotic velocity lift. 

The other approach we shall consider here is known as the modified Newtonian dynamics (MOND), 
which was proposed by Milgrom (1983); in this case there is no dark matter at all, rather a deviation from 
the usual Newton's second law would occur when one is dealing with very 'small' accelerations, where 'small' 
means small with respect to some (in principle, universal) critical acceleration parameter, ao. Instead of 
F = m g a, one would write: 



g N = n(a/a )a, (Bl) 

where gjv is the conventional gravitational acceleration, a is the true acceleration of a particle with respect 
to some fundamental frame (a = |a|), and /i is a function of a/ao of which only the asymptotic forms 
fi(a/ao 3> 1) ~ 1, fi(a/ao <C 1) ~ a/ao are known. It can then be seen that for accelerations much larger 
than the acceleration parameter ao, fJ> « 1 and we recover the Newtonian dynamics. For the rotation law, 
the usual expression remains to be valid: v^/r = a, and also gN ~ MGr~ 2 ; combining this with equation 



(Bl), we get the asymptotic velocity: 



v% = GMa . 



In this case the acceleration parameter ao can be taken as fixed so there is only one free parameter, again 
the ratio M/L. Alternatively, ao can also be considered as a free parameter. Although this approach works 
very well when fitting the rotation curves, there is no a priori reason to believe that a deviation of this kind 
could indeed occur, so this is usually considered to be a purely phenomenological description. 
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